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Two different approaches to the Roberts dynamo problem are considered. Firstly, the 
equations governing the magnetic field are specified to both harmonic and subharmonic 
solutions and reduced to matrix eigenvalue problems, which are solved numerically. Sec- 
ondly, a mean magnetic field is defined by averaging over proper areas, corresponding 
equations are derived, in which the induction effect of the flow occurs essentially as an 
anisotropic alpha-effect, and they are solved analytically. In order to check the reliability 
of the statements on the Karlsruhe experiment which have been made on the basis of a 
mean-field theory, analogous statements are derived for a rectangular dynamo box con- 
taining 50 Roberts cells, and they are compared with the direct solutions of the eigenvalue 
problem mentioned. Some shortcomings of the simple mean-field theory are revealed. 



1. Introduction. The dynamo model proposed by Roberts 1972 [13] has 
been chosen as the starting point for an experimental demonstration of homoge- 
neous fluid dynamo at the Forschungszentrum Karlsruhe [H [3] 21 [6l [TO] [TTJ [17l [HI 
[213] . The fluid velocity field considered by Roberts which is of particular interest 
in this context is given by 

U = U(sm(Y/Lu), cos(X/Lu), x(cos(X/Lu) - cos(Y/Lu))) ■ (1) 

Here a Cartesian coordinate system (X, Y, Z) is used. The flow pattern is sketched 
in Fig. [T] Ljj is the length of the diagonal of a cell in the AY-plane and the 
parameter \i which is a constant, determines the Z-component of the flow and 
so the helicity of the velocity field. Roberts has demonstrated that a flow of this 
kind is capable of dynamo action. He investigated, however, only magnetic fields 
which show the same periodicity in X and Y as the flow pattern. These fields, 
which we call here "harmonic fields" , possess parts which do not depend on X 
and Y but only on Z or, in other words, they have infinite wave lengths in the 
X and Y directions. As Roberts himself pointed out the considered flow allows 
also non-decaying magnetic fields with finite wave lengths in all directions. For 
a particular case such fields were investigated by Tilgner and Busse [19 , who 
called them "subharmonic" , and in a more general frame by Plunian and Radler 
[5]. Despite the finite dimensions of the Karlsruhe experimental device many 
estimates concerning excitation conditions etc. have been made on the basis of 
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Figure 1: The stream lines of the Roberts flow in the zy-plane. They coincide with 
the isolines of the velocity in z-direction. Dimcnsionless coordinates x = X/Ljj 
and y = Y/Ljj are used. 



findings about harmonic magnetic fields. It is, however, of high interest to compare 
these results with such derived from results on subharmonic fields. 

In this paper we start with the basic equations of the Roberts dynamo problem 
and some general consequences (Section [5]), present some findings on its harmonic 
solutions (Section[3]) and explain a mean-field approach to the dynamo problem on 
that level (Section [4j. After that we turn to subharmonic solutions and give some 
results for them (Scction[5|). We then deal with the Karlsruhe experiment, derive in 
the framework of a mean-field approach and under simplifying assumptions on the 
dynamo module an excitation condition and compare it with a corresponding result 
of the subharmonic analysis (Section[5]). Finally summarize the main consequences 
of our findings (Section [7]). 



2. The Roberts dynamo. To discuss the Roberts dynamo problem in 
some detail we consider the induction equation governing the magnetic field B, 
assuming that it applies in all infinite space. We use its dimcnsionless form 

— = V x (u x B) + RmT^B, V • B = (2) 

with 

u = (smy, sinx, x(cos x — cosy)) . (3) 

Instead of X = (X, Y, Z) we have introduced here dimensionless coordinates x = 
(x,y,z) defined by x — X/Ljj, instead of U the dimensionless velocity u defined 
by u — U/U, and we measure the time in units of Ltj/U. Further R m is the 
magnetic Reynolds number defined by 

Rm = ^ (4) 

V 

with 7] being the magnetic diffusivity of the fluid. 
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For a steady flow as envisaged here we may expect solutions B varying like 
exp(pi) in time, where the real part of p is the dimensionless growth rate. In 
this case an eigenvalue problem for B with the eigenvalue parameter p occurs. 
Furthermore, since the flow is z-independent, B can be assumed to possess the 
form 

B = R{b(x, y, k) exp(pt + ikz)} , (5) 

where b(x,y, k) is a complex vector field independent of z, and k a dimensionless 
wave number with respect to the z-direction. When inserting ([5]) into @ we find 

pb + {u ■ V)6 = (b ■ V)« - \ku z b + i? m _1 (V 2 - k 2 )b (6) 

V • b + ikb z = 0. (7) 

The x and y components of © are equations for b x and b y which do not contain 
b z . They constitute the mentioned eigenvalue problem. After solving it, we can 
calculate b z from ([7]) without any integration. 

We may easily conclude from ^ that the results for p for any x can be 
inferred from those for \ = 1 with the help of the relation 

p{k, R m , X ) = p(k X , Rm, 1) + k 2 R m -\ X 2 - 1). (8) 

As long as we deal with direct solutions of the Roberts dynamo problem (up to 
Section [5]) we therefore restrict our attention to the solutions with X = 1. In the 
discussion of the Karlsruhe dynamo experiment (Section [6]) we admit also other 
values of x- 

3. Harmonic solutions. As mentioned above, Roberts solved the rele- 
vant equations only for magnetic fields with the same periods in x and y as the 
flow pattern. We consider first this case only, in which we speak of "harmonic 
solutions". Then b{x,y, k) must have the same periodicity in x and y as the flow 
pattern. 

Solving the eigenvalue problem defined by © Roberts found growth rates, 
that is real parts of p, in its dependence on k for values of R m up to 64 as shown 
in Fig. [5] The imaginary parts of p proved to be equal to zero (numerically 
always close to zero), attesting that the dynamo instability is an absolute one. 
This implies that the magnetic field geometry is stationary while the intensity in 
general varies in time. 

Soward |14j has shown that in the limit of large R m the order of the maximum 
of the dimensionless growth rate, p max , is given by p max = (9(ln(lni? m )/lni? m ). 
It occurs at a wave number & max for which fc max = 0(R m 1 ^ 2 (\nR m )^ 1 ^ 2 ). That 
is, Pmax — * as R m — * oo. Thus the Roberts dynamo proves to be a slow one. 
This applies not only in the context of harmonic solutions, for it turns out that 
other solutions never grow faster than the fastest of the harmonic ones. 

Results for finite R m obtained in our calculations are given in Fig. [2] and in 
Tabled] see also [HI [16]. Note that the quantities Ti = p max lni? m /ln(lni? m ) 
and T2 = fc max (lni? m ) 1 / 2 /i? m 1 ^ 2 given in Table [1] approach constant values as R m 
grows and thus illustrate the mentioned asymptotic laws. 
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Figure 2: The growth rate p versus k for various R m (given by the labels at the 
curves) and \ = 1. For each given R m the circle corresponds to the maximum 
growth rate p max occurring at fc max . 
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Table 1: Maximum growth rates p m ax, the corresponding wave numbers k m ax and 
the quantities T\ and T2 for various R m and x = 1. 
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4. A mean— field approach. The Roberts dynamo with harmonic mag- 
netic fields B can also be understood in terms of mean fields. For any field Q we 
define a mean field Q by averaging over an area of one periodic unit, or four half 
cells, of the flow pattern in the xy-plane as depicted in Fig. [TJ for example 

/ + 7T />+7T 
/ Q(x + t,y + v, z ,t) d t d v- (9) 
-7T J — 7T 

Whereas B is in general non-zero we have u = 0. 

Subjecting the induction equation @ to this kind of averaging we obtain 

— = V x £ + Rrn-^B, V-B = 0, (10) 
at 

where £ is a mean electromotive force due to the fluid motion, 



e = uxB. (ii) 

We admit here no other magnetic fields B than harmonic ones in the above sense. 
Then B and £ depend no longer on x and y but only on z and t. 

Let us further use Fourier representations for B, B and £ according to 



Q(x,y,z,t) = J Q(x,y,k,t)e l " z dk, (12) 

where the integral is over — oo < k < oo. The corresponding representation of 
B clearly includes the ansatz ([5|). B depends on x, y, k and t, but B and £ 
depend only on k and t. The requirements that B, B and £ are real lead to 

B (x, y, k, t) — B(x, y, — k, i) and analogous relations for B and £. 
According to (JTTJ) we have 



£{k,t) = u{x,y)x B{x,y,k,t). (13) 

With standard reasoning of mean-field theory we conclude that £ is linear and 
homogeneous in B. For the sake of simplicity we assume that £ at a given time 
depends only on B at the same time, that is, ignore any dependence on B at 
earlier times. We therefore write 

£ i {k,t) = a ij {k)B~ j (k,i), (14) 

where &ij is a complex tensor determined by the fluid flow. Analogous to £ and 
B it has to satisfy a*j(k) = a.ij{— k). From the symmetry properties of the u-field 
we conclude further that the relation (I14p remains its validity if both £ and B 
are simultaneously subject to a 90° rotation about the z-axis, and that £ has no 
z-component and does not depend on the z-component of B. The first fact means 
that the tensor dy is axisymmetric with respect to the z-axis, that is, can only 
be a linear combination of 5y, aej and etjiei, where 5y is the Kronecker tensor, 
6iji the Levi-Civita tensor and e the unit vector in z-direction. The second fact 
implies that Sij and eiej occur only in the combination 6y — eiej. Thus we may 
write 

&ij(k) = -ai(k)(Sij - eiCj) + ika 2 {k)e i:j iei (15) 
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with two complex functions di(fc) and &2(k). The signs in this relation and the 
factor ik of the last term have here to be considered as arbitrary but will prove to 
be useful in the following. 

As can be easily seen from the equations (fTU)) and (fTTj) their solutions must 
have the form 

B = $t{Cexp(pt + ikz)} (16) 

with a complex vector C lying in the xy-plane. If we use (|14[) and (|15[) we find 
that 

C x = ±iC y (17) 

and 

p + R m ~ x k 2 = a{k)k, a(k) = ±&i(k) - a 2 (k)k (18) 

The two signs in (jTTJ) and (fl~8|) indicate the existence of two classes of solutions 
and had been already identified by Roberts [T3] ■ Relations of the same kind have 
also been found by Soward 14] (who used the notations ac and old instead of -6t\ 
and 0.2)- We can calculate a(k) with the help of the harmonic solutions of the 
induction equation @ obtained with the ansatz ([5]). There are again two classes of 
these harmonic solutions, which correspond to the two signs in (fl~7|) and (fT8| . The 
two classes lead to two different functions a(k). In agreement with the finding 
that p is real also a proves to be real. Among the two function a(k) only the 
largest one (the upper sign in (|17[) and (|18p) has been retained as it corresponds 
to the largest growthrate p (see also [5])- Fig. [3] shows the dependence of a on k 
for various R m and % = 1. The values of a for arbitrary \ can be inferred from 
the values for \ — 1 with the help of the relation 

a(k,R m ,x) = Xa( k X,Rm,l) , (19) 

which follows from (jHJl and (fT5|) . 

Since a is real we may assume that also 6l\ and 012 are real. Moreover we 
may then conclude from &*j(k) = (— fc) that both are even in k, that is a\{k) — 
&i(—k) and &2(k) = a2(—k). If we know d(fc) for positive and negative k we may 
find di(fc) and d2(fe). 

We may conclude with the help of (Tl4]) and (fTS]) that 

= - J &i(k)(B~(k,t)-(e-B~(k,t))e)exTp(ikz)dk 
d f 

-ex — a 2 (k)B(k,t) exp(ifcz)dfc . (20) 

This is equivalent to 

5(;M) = -^y«i(C)(S(z + C,<)-(e-B(z + C,i))e)dC 

1 a 

e x — 

2tt 9z 



-T^ex-?- I a2(()B(z + C,t)d( (21 i 



with 



ai(C) = y di(fc)exp(ifcC)dfc, a 2 (C) = / d 2 (fc) exp(ifcC)dfc . (22) 
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Figure 3: The coefficient a versus k for various R m (given by the labels on the 
curves) and x = 1- 



The integrals are again over — oo < k < oo or over — oo < ( < oo. Note that both 
a.1 and a-i are even in £. 

Let us now expand &ij(k) in a Taylor series. ^From (|15p we have 



<%(&) = -ai(0)(5. iJ - CiCj) + i/ca 2 (0)e y 7ez + 



(23) 

The second relation (|18|) together with the symmetry properties of ai and d2 leads 
to 

a x (0)=a(0), a 2 (0) = -^(0). (24) 



The corresponding expansion of £ reads 



dB 



-a±(B - (e • B)e) - 0e x — + 



where 



a± 



d(o) 



C»2 



(25) 



(26) 



The same result can be derived from (|2T|) by expanding B(z + £, i) with respect to 
£. The first term on the right-hand side of (|23|) describes the anisotropic a-effect. 

Since here e x = V x B it seems at the first glance reasonable to consider the 
second term as a contribution to a mean-field diffusivity. By a reason mentioned 
later in Section [6.21 however, this interpretation is not compelling. 

Values of a ± and (3 for various R m are given in Table [5J Remarkably enough 
(3 is negative as soon as R m exceeds a value between 1 and 2. This means that 
then the corresponding term in (|25p supports the dynamo action of the flow. For 
the limit of large R m it was shown that a± = 0(R m ~ 1 ^ 2 ) [2H]. This is illustrated 
by the values of the quantity T 3 = a^R„ } l l' 1 given in Tabled 



7 





0.25 


0.5 


1 


2 


4 


8 


16 


32 


64 


a ± 


0.24 


0.44 


0.66 


0.62 


0.39 


0.27 


0.19 


0.13 


0.09 


P 


0.25 


0.46 


0.54 


-0.3 


-0.4 


-0.12 


-0.07 


-0.05 


-0.03 


r 3 


0.12 


0.32 


0.66 


0.87 


0.79 


0.75 


0.75 


0.75 


0.75 



Table 2: The coefficients a± and (3 of the expansion (f2"S"]) and the quantity 1^ for 
various R m and \ = 1. 

5. Subharmonic solutions. In order to check a simple mean-field the- 
ory of the Karlsruhe dynamo experiment we are interested in solutions B of the 
induction equation ^ with period lengths exceeding those of the flow pattern, 
which we call "subharmonic solutions" . We consider the case in which the period 
lengths of B are larger by an integer factor N than those of the flow pattern. As 
mentioned above, this problem has already been investigated by Tilgner & Bussc 
[I!?] for a few special values of ./V and later by Plunian and Radler 5 . 

We focus our attention again on the induction equation ^ governing the 
magnetic field B in all space. We use again ([5]) but consider b no longer as a 
field with the same periodicity in x and y as the flow pattern. Instead we put 
b = b(x, y, f x , f y , k) exp(i(f x x + f y y)), where b has now the same periodicity as 
the flow pattern and f x and f y are subharmonic wave numbers in the x and y 
directions. In that sense we look for solutions of the induction equation of the 
form 

B = di (b(x, y, f) exp(pt + if ■ xfj . (27) 

with b(x,y, f) being a complex periodic vector field with the same period length 
in x and y direction as the flow pattern, / the real vector (f x , f y ,k), and p again 
a complex quantity; for more details see [5]. We restrict our attention here to the 
case fx = fy — f= l/N. Then the period lengths of the magnetic field B are 
just N times that of the flow pattern. The harmonic solutions B discussed above 
correspond to the limit N — > oo. 
Inserting (|27| in ([2]) we find 

pb + {u ■ V)6 = (6 • V)tt - i(f ■ u)b + R^ 1 (V 2 b + 2i{f ■ V)6 - fb), (28) 

V b + ib- f = 0. (29) 

The system (|28j) defines an eigenvalue problem with p being the eigenvalue pa- 
rameter^- It has been solved numerically. Marginal values of R m versus k for 
fx = fy = / are shown in Fig. [4] for different values of / and again x = 1. For 
/ 7^ there are both a critical R m and a critical k below which dynamo action is 
not possible. In general p is no longer real, that is, we have no longer stationary but 
moving field structures. We point out that relation jS]) allows us the calculation 
of p for arbitrary x from that for \ = 1. 

6. Applications to the Karlsruhe experiment. 

lr These equations have already been derived by Roberts (1972) though used only with f x = 
f y = in his numerical calculations. 
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Figure 4: Marginal values of R m versus k for f x = /„ = / and x = 1. The curve 
(a) belongs to / = 0, (b) to / = 1/10, (c) to / = 1/8, (d) to / = 1/6, (e) to 
/ = 1/4. 

6.1. The Karlsruhe "dynamo module" . The essential piece of the Karlsruhe 
dynamo experiment is the "dynamo module" , a cylindrical container with both 
radius and height somewhat less than 1 m, through which liquid sodium is driven 
by external pumps [3] 01 \17\ [18] . By means of a system of channels with conducting 
walls, constituting 52 "spin generators", a helical motion is organized. The flow 
pattern is similar to that defined by ([1]). The 52 spin generators correspond to 26 
periodic units of the flow pattern such as the one shown in Fig. [1] The arrangement 
of the pumps allows to vary the parameters U, or R m , and x independently from 
each other. 

6.2. A simple mean-held theory of the experiment. In order to give an es- 
timate for the self-excitation condition of the experimental device a simple mean- 
field theory has been developed [7l [H [9[ [10l [IT] [12]. The mean magnetic field 
B defined as above is assumed to satisfy the equations (|T0|) and (fTTj) inside the 
dynamo module and to continue in some way in outer space. Of course, in this 
context B can no longer be independent on x and y, and £ can no longer have 
the simple forms (|20l) . (f2"Tj) or (f2"5|) . Relying on some traditional concept it was 
assumed that the variations of B in space are sufficiently weak so that £ in a 
given point can be represented by B and its first spatial derivatives in this point. 
Together with the symmetry properties of the Roberts flow this leads to 



where ojj_, /3_l> P\\ and /?3 are constants depending on R m and x> an d e is again 
the unit vector in z-direction JS[ [TTJ . As in (f2"S"l) the term with a± describes the 



£ = 



a±(B - (e ■ B)e) - /3j_V x B — — /3±)(e ■ (V x B))e 
3 ex (V(e.B) + (c-V)B), 



(30) 
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anisotropic a-effect acting in the xy-plane only. The terms with f3± and /3ii can 
be interpreted by introducing an anisotropic mean-field diffusivity different from 
the molecular magnetic diffusivity. Finally the term with /?3 describes a part of £ 
depending on derivatives of B which cannot be expressed by V x B and therefore 
not be interpreted as a contribution to a modified diffusivity. Several results have 
been derived on the dependence of a± on the fluid flow [51 [7J [HI M El HI] , an d 
also such on (3±, /3|| and /?3 [6]. 

For a field B not depending on x and y and having no z-component we have 
V x B = e x and the last three terms on the right-hand side of (|30j) can be 

written in the form — (/3± + /3z)e x As to be expected, in this special case the 
structures of ([25)) and ([30)) coincide, and we have (3 = /3± + (3 3 . Our above remark 
on the /33-term in ([30)) explains why the interpretation of the /3-term in ([25]) as a 
contribution to a mean-field diffusivity is not compelling. 

The assumption on small variations of B in space means in particular that 
B does not change markedly across a spin generator. In that sense the usage 
of (|30p in a theory of the dynamo module can only be justified for a very large 
number of spin generators within the module. Quite a few solutions of the equation 
(fTU)) for B, applied to the dynamo module, with £ according to ([30]) and various 
boundary conditions have been calculated [HI UHl [XX] ■ In most cases, however, no 
other contribution to £ than the a-effect, that is, only the first term on the right- 
hand side of pop was taken into account. Contributions with higher than first 
derivatives of B have never been considered. By these and other reasons a check 
of the results of the simple mean-field theory on a way that avoids the mentioned 
shortcomings seems very desirable. 

6.3. Comparison of results of mean-field approach and subharmonic analysis. 
For this purpose we deal now with a very simple model of the dynamo module. 
We consider no longer a cylindrical module but instead a rectangular "dynamo 
box" with a quadratic basis area in the xy-pl&ne and denote the edge lengths 
of the box in this plane by L and its hight by H . Thinking of the shape of the 
real dynamo module we put L/H = 2. We will study the excitation condition 
for a mean magnetic field B which satisfies the equation (fTU)) and the relation 
(jlip in all space and is periodic in x and y with the period length 2L and in z 
with the period length 2H . This periodicity means that the dynamo box contains 
just a "half wave" of the field B. For the sake of simplicity we use ([TT]) in its 
reduced form containing no other induction effect than the a-effect. We will 
then compare this excitation condition with that for a subharmonic _B-field whose 
longest wave lengths show the same periodicity, that is, which fits in the same 
sense to the dynamo box. In this context we put N = 10 so that an area of 
(2L) 2 in the xy-plane contains just 100 period units, that is, 200 cells of the flow 
pattern, consequently the basis area of the dynamo box 50 cells, which have to be 
compared with the 52 spin generators in the real dynamo module. N = 10 means 
/ = 1/10, and with L/H = 2 we arrive at k = 1/5. 

Instead of a realistic boundary condition for the dynamo module we use here in 
fact the condition of periodic continuation of the magnetic fields both on the mean- 
field and the subharmonic level. Such a condition might be in general problematic 
but seems acceptable for the comparison which we have in mind. 
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From Fig. [3] we see that the value of a for k = 1/5 can, except for small 
R m , not be inferred from a(0) and §f (0) only. This suggests that there will be 
discrepancies between the excitation conditions obtained with a mean-field theory 
which ignores contributions to S with higher than first-order spatial derivatives 
and those derived from the subharmonic analysis. 

As already explained we assume for our mean-field consideration that equa- 
tion (p~0|) and the reduced form of ([30)) . that is, 

d~B 

dt 



= -Vx(a i (B-(e.B)e))lt 1 V 2 B, V-B = 0, (31) 



apply in all space with constant a±. We may represent B as a sum of a poloidal 
and a toroidal part, 

B = -Vx(exVS)-exVT, (32) 

with two defining scalars S and T. Inserting this in (|3ip and dropping unimportant 
constants we find 

R m - 1 AS-a±T- ^ = 
at 

R m ^AT + a ± —-— = 0. (33) 

The special periodic solution of (f3"Tj) which we are looking for is obtained with the 
ansatz 

S = So cos (fx) cos (fy) cos (kz) exp(pi) 

T = (T /S )S, (34) 

where So and To ar e constants, / and k the parameters specified above and p, 
which will prove to be real, is again the growth rate. When inserting this in ([33]) 
we arrive at two linear homogeneous equations for So and To. The requirement 
that they allow non-trivial solutions leads to 

P=- V ±\a ± \k. (35) 

-Km 

Growing B are possible in the case of the upper sign of the last term if is 
sufficiently large. The excitation condition reads 

\«±\>^rf- (36) 

In the representations of results on a± on which we now rely [SJ [71 [HI [S] 
the latter is given in its original dimension so that it corresponds U a± with our 
dimensionless a± . Furthermore these results are given in terms of the two magnetic 
Reynolds numbers i? m i| and R m ± for the flow in the xy-plane and in z-direction, 
respectively. These are connected with our R m and \ by 

Rm\\ = ——XRm, Rrn± = 2R m . (37) 
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The marginal states of the dynamo, in which B neither grows nor decays, 
are given by pairs of R m ± and R m \\, or by the corresponding neutral curve in the 
-Rm±-Rm|| - diagram. We may represent the result for arbitrary / and k by using 
the modified magnetic Reynolds number R m *\ defined by 

Rm ^ \^M Rm ^W+k^ (38) 

instead of R m \\- Note that R m *\ is no longer determined by U alone but also by 
B. Fig. [5] shows a i? m j_i? m | -diagram in which curve (a) gives just the result of 
our mean-field calculation. Clearly dynamo action requires that i? m |j exceeds a 
critical value, ft appears, however, to be possible for any R m ± if only Rmu is 
sufficiently large. 

Let us now compare this result with that for a corresponding subharmonic 
solution B of the induction equation. In Fig. [5] the curve (b) is the neutral one 
for the subharmonic solution with the values of / and k specified above. Clearly 
dynamo action requires now not only that R m *\ exceeds a critical value but also 
that R m ± lies above such a value. In addition for each given R m ± allowing dynamo 
action the marginal value of i? m jj derived in the subharmonic approach is higher 
than that concluded from the mean-field approach. In the range of R m ± between 
1.2 and 2, which corresponds to the actual situation in the Karlsruhe experiment, 
the deviation is larger than 20 %. Of course it will become smaller in a comparison 
with a mean-field model which involves also the induction effects connected with 
first derivatives of B indicated in (|30|) : see [6]. But even then the mean-field 
approach underestimates the requirements for self-excitation. 

7. Conclusions. We have dealt with several aspects of the Roberts dy- 
namo problem and derived some results which are of interest for the Karlsruhe dy- 
namo experiment. Although a rectangular dynamo box was considered, there are 
good reasons to assume that the main conclusions apply as well to the real exper- 
imental device with a cylindrical dynamo module. In the framework of the simple 
mean-field theory of the experiment self-excitation seems possible for arbitrary 
values of the magnetic Reynolds number R m ± describing the flow perpendicular 
to the axes of the spin-generators if only the magnetic Reynolds number R m || for 
the axial flow is sufficiently large. An analysis based on subharmonic solutions re- 
vealed that a dynamo is only possible if both R mJ _ and R m ii exceed critical values. 
Apart from this it was found that the simple mean-field theory underestimates 
the excitation condition of the dynamo. This discrepancy of the mean-field results 
with those obtained with subharmonic solutions cannot be completely removed by 
taking into account the effect of the mean-field diffusivity. 

The authors are indebted to Dr. M. Rheinhardt for several helpful comments 
on the draft of this paper. 
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